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We present full numerical solutions to the system of a global string embedded in a six-dimensional 



space time. The solutions are regular everywhere and do confine gravity in our four-dimensional 
world. They depend on the value of the (negative) cosmological constant in the bulk and on the 
parameters of the Higgs potential, and we perform a systematic study to determine their allowed 
values. We also comment on the relation of our results with previous studies on the same subject 
and on their phenomenological viability. 
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I. INTRODUCTION 



Field theory constructions in which our four-dimensional world consists of a topological 
defect embedded in a higher dimensional space-time are more and more frequent these days, 



m. 



Such defect could either 



although they were originally proposed twenty years ago 
be a domain wall (if the total number of dimensions is five), a string (six), a monopole 
(seven) or an instanton (eight). 

The gravitational field of these topological defects becomes relevant in this context. In 
D = 4 that of the domain wall has been thoroughly studied . Whereas it 

seems to be non-static from the perspective of an observer on the wall, it is actually the 
domain wall which is non-static, moving in a static Minkowski space-time. The gravitational 
field of the global monopole is as well static and well-defined [lOj], and it is only the global 
string which happens to give rise to a static and singular metric outside the core of the defect 
jrH . It was later realised that this singularity could be cured by adding time dependence to 
the metric jl2|. 

Coming back to higher dimensions, it is then mandatory to explore the gravitational field 
generated by these defects if one is willing to build up realistic models in which our four 
dimensional world is to become one of them. Following the work of Randall and Sundrum 
jl^ . the domain wall in D = 5 has been extensively studied, and we have nothing new 
to add to it. Also, the generalization of the global monopole to higher dimensions, which 
results in a D = 7 space-time, has been studied jl^j], resulting in, again, a well-defined static 
metric. The issue then remained of whether the gravitational field generated by the global 
string would still result singular in its D = 6 format. 

This issue was first investigated by Cohen and Kaplan jisj . who concluded that, similarly 
to the four- dimensional case, the metric around the global string is still singular. Then 
Gregory argued that, again analogously to the D = 4 case, adding time dependence to 
the metric should cure that problem and that an equivalent procedure would be to add to the 
static metric a negative cosmological constant in the bulk. Higher dimensional extensions of 
these claims were presented in |l7[. In these last two articles analytic arguments were given 
to support them, although a full numerical solution was not presented. 

This is precisely the main aim of the present letter. We want to present numerical 
evidence of the existence of solutions that confine gravity around a global string in which 
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core our four-dimensional world exists. The geometry of the D = 6 space-time is that of a 
cigar, and it behaves asymptotically as AdS§ x Si. These solutions are numerically hard to 
obtain and rely on a very precise tuning of the value of the (negative) cosmological constant 
in the bulk. 

In the next section we describe the model we are going to work with, namely the system of 
a scalar field together with gravity in D = 6, and we write the set of equations to which their 
dynamics reduce. In section 3 we show our solutions, explaining the technical involvement 
of finding them, and also their physical interpretation. We compare with related models 
already published in the literature and, in section 4, we conclude. 

II. THE MODEL 

A global string is a topologically non trivial solution of the following Lagrangian 

C = ^{V A $)tV A $-V($) (1) 

where the potential has to exhibit a global U(l) symmetry. In particular, we shall study the 
so-called Mexican-hat potential, i.e. 

= A($ t $ - v 2 ) 2 , (2) 

Our goal is to find numerical solutions for the equations that govern the system formed by 
the field when considered in the context of the six-dimensional geometry defined by 

ds 2 = M 2 {r)7 1 ^x u - dr 2 - L 2 {r)d9 2 , (3) 

where (r, 9) are the coordinates of the transverse space and M(r), L(r) are the warp factors. 
In particular we are looking for solutions that would result in gravity trapping in four space- 
time dimensions. 

Our starting point will be the action for the D = 6 system 

s =-/ tf ^S + A " £ ) • (4) 

with g the determinant of the metric, eq. Q, k 2 = 1/Mq, and M 6 is the D = 6 Planck mass. 
A is the bulk cosmological constant and C is given by eq. ((H). 
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The Einstein equations for this system are given by 

1 1 

Rab ~ ^9abR = J^(^9ab + T ah ) , (5) 

where a = 1, ... ,6 and the equation of motion for the field $ is given by 

BV 

9 ab VaV b $ + w = . (6) 

In order to simplify the presentation of results we shall parametrize the scalar field as 

$ = vF{r)e iqd , (7) 

and we define m(r) = M'{r)/M{r) = dha.M(r)/dr. Our task is, therefore, to solve a set of 
second order differential equations for the variables f(r), L(r) and m(r) for certain values 
of the parameters A, v and q. The actual equations, already assuming Mq = 1, are given by 



J a Ji /„. 2zt>/2 „2„,2p2 \ 

±- + 3m' + 6m 2 + 3^m = -A - ( + + \ V \ F 2 - l) 2 ) 

L L \ I 1L Z J 



v 2 F' 2 


q 2 v 2 F 2 


2 


+ 

2L 2 


v 2 F' 2 


q 2 v 2 F 2 


2 


2L 2 



Am' + 10m 2 = -A - ( + \v\F 2 - l) 2 J , (8) 

2 ^ m + 3m 2 = -A - - (-v 2 F' 2 + + \v\F 2 - l) 2 



L 2 V 2L 2 

F" + ( 4m + ^ V' = 4AFt; 2 (F 2 - 1) 1 ^ 



L / v L 2 ' 

Therefore we have to solve three differential equations (the third one of the previous system 
is a constraint) with a mixture of boundary conditions defined at the origin and at infinity. 
This is a typical boundary value problem. 

To be more precise, let us elaborate on the choice of boundary conditions: at the origin 
(r — > 0), we demand a regular geometry at the core of the string and the absence of deficit 
angle in our solution. This translates into 

m(0) = 0, L(0) = ,L'(0) = 1 . (9) 

Moreover, a local analysis by power series shows that, near the origin, 

F{r) = , 

L(r) = r + l 3 r 3 , (10) 
m{r) = mir . 
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Substituting these ansatze in the previous equations (JBJ) we get 

mi = -\(A + V(0)) , 
h = ^(A + V(0)-2f 1 ) . (11) 

fx remains a free (shooting) parameter. 

Far away from the core, at r — * oo, we demand that all three functions in eqs. go to 
constants. The metric in this region is then assumed to be cigar-like 

dsl, = e 2m °° r r] lu ,x^x v - dr 2 - L^dO 2 . (12) 
Substituting again in eqs. (jSJ), and reinstating factors of M 6 where necessary, we get 

= ^^t^ . (13) 

For phenomenological reasons, i.e. in order to have gravity trapping in D = 4, we are looking 
for solutions that correspond to < 0. The second warp factor is given by 



°° 2M||m 00 | 



(14) 



whereas is obtained by solving the equation 

A + V(f oo )-2f oo V'(f oo ) = . (15) 

The first interesting conclusion to be drawn is that the solution of the previous equation will 
not correspond to the field settling at the minimum of its potential, i.e. V'(foo) — 0. The 
presence of the negative cosmological constant, in other words, the interplay of the scalar 
field with gravity induces the field to settle just before reaching its minimum, as we will see 
next. 



III. RESULTS 

Once we have defined the system we want to work with, we can look for solutions that 
satisfy our requirements. In order to perform a numerical analysis we have used a relaxation 
method which would look for a solution once values for A, q and v were specified. This 
means that we replace the system given by eq. (JEJ) by a set of finite-difference equations 
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(FDEs) on a mesh of points that spans from zero to a sufficiently large value for r. Given 
that we are starting off with a system of 5 coupled first-order equations (remember that, in 
eq. (JHJ), one equation is a constraint and another is already first order) represented by FDEs 
on a mesh of M points, the solution consists of values for 5 x M variables. The relaxation 
method determines the solution by starting with a guess and improving it iteratively. As 
the iterations improve the solution, the result is said to relax to the true one. 

From now on we will focus on the case q — 1. Then, for every value of v there is a 
unique value of A (which we shall call A c ) which gives us a regular (with no deficit of solid 
angle) solution everywhere. This was already pointed out by Ruth Gregory and here 
we present the first numerical evidence of this statement. 

The way in which we determine A c is as follows: in our relaxation code we give as 
boundary conditions the values of F(r), L(r) and m{r) at the origin, and those of F'{r) 
and L'(r) at infinity, which corresponds to a typical boundary value problem. This means 
that we do not know a priori the value of L'(0), which determines whether or not there is 
a deficit of solid angle. Therefore, for every value of v we try different values of A until we 
find a solution with L'(0) = 1. 

In order to simplify the description of the results, we shall use adimensional variables to 
draw the plots, which requires the following reparametrization. 



This already shows that the results are independent of the value of A. 

In figure 1 we show the evolution of the scalar field, represented by f(x) normalized to its 
value at infinity, which is given by f^, = 0.99577, i.e. the field gets stabilized just before its 
minimum located at v — 1. As for the two warp factors, L(x) and m(x), shown in figures 2 
and 3 respectively, we can see how both go to constants at infinity. The associated geometry 
of the space-time is, therefore, that of a cigar, which goes to AdS^ x S\ asymptotically. 

We have performed a systematic analysis of the parameter space, by varying the value of 
v and adjusting the corresponding one of A c . The results are shown in figure 4, where we 
plot the warp factor m(x) as a function of x for v = 1.4, 1, 0.7, 0.6. As we can see, the smaller 



x 





(16) 
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FIG. 1: Plot of /0)//oc as a function of x for v = 1, A c = -0.0671. 
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FIG. 2: Plot of L(x)/Loo as a function of x for u = 1, A c = -0.0671. 

v is the longer it takes the gravity fields to settle at their asymptotic values. This is just 
reflecting the fact that, for smaller v, the scalar field tends to settle closer and closer to its 
minimum, playing less of a role in the stabilization of the warp factors. In other words the 
problem becomes a typical two-scale one, where the scalar field quickly converges whereas 
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FIG. 3: Plot of m(x)/m OQ as a function of x for v = 1, A c = —0.0671. 




FIG. 4: Plot of m(x)/m OQ as a function of x for v = 1.4, v = 1, v = 0.7 and ?> = 0.6 



To further comment on numerical issues, let us describe how difficult it is to obtain these 
solutions. As mentioned before, for a given model, there exists a regular solution just for 
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a unique value of A. This was thoroughly explained in Ref. jltsj . and here we shall merely 
repeat the main arguments and give a numerical proof. Essentially it is fair to study this 
system in its asymptotic region by assuming the scalar field f(r) to be at its minimum and 
working with just the two warp factors L(r) and m(r). Then we can define the following 
variables 



1 / V 

x = 3m. H 

-2A I L 



1 V 

y 



(17) 



-2AL ' 

which, together with the independent variable p = x\/—2A, define the following autonomous 
dynamical system, 



xy 4 2 

(18) 



x = y 

3 3 y 



, 4 2 3 

V = ^x{x -y)-y - - . 

The primes here mean derivatives with respect to p. 

In order to understand the structure of these numerical solutions, we must calculate the 
critical points of this system. Those are given by 

cf = (±11' °) saddle 

(19) 




= 5 ' y/4 5 ' 

The solutions we have found, which we have shown in previous plots, correspond to flowing 
towards c[ and it is now easy to understand why they are so hard to obtain: this is a saddle 
point, which is next to a repellor, given by (which, by the way, would be the critical 
point describing a geometry of the type AdSe). Essentially only one trajectory, the one 
corresponding to A c , ends up in c~[ which can be matched to a regular solution near the core 
of the string. Solutions ending up in c\ would correspond to a four dimensional metric that 
would blow-up. 

For completeness we give a plot of the phase space for this autonomous system in figure 5, 
where we have also inserted the solution we have found (with a thick line). As we can see, 
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FIG. 5: Plot of the (x,y) phase space, for v = 1. 



far enough from the saddle point a[, the scalar field F(r) still plays a role in the evolution 
of x and y and, therefore, the inserted trajectory does not fit too well within those obtained 
from solving the system (|18|). We have highlighted this by using dashed lines in that part 
of our solution. As the field F(r) plays less of a role, the trajectory starts to fit 'naturally' 
into the phase space drawn, and that is represented by the continuous thick line. 

In fact, one can even go further and obtain an analytic expression for the two warp 
factors L(r) and m(r) as they approach their asymptotic limit. This is done by rewriting 
the autonomous system Eq. ()18j1 in the linear approximation around c~[ , and solving for the 
two functions. The final result is that 



The A factor, which determines the normalization of the curves, can then be extracted from 
our numerical results. We have checked that the results we get for A from the two curves, 
i.e. L(r) and m(r), are compatible with each other and that the fit to both functions in the 
asymptotic region is extremely accurate. 

Next, we would like to discuss the role played by the cosmological constant in the bulk, 




(20) 
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A. As it was mentioned at the beginning of this section, we have numerically checked that, 
for every value of v there is a unique value of A, which we denote as A c , that gives us a 
regular solution everywhere. We have explored the parameter space defined by A c / A and v, 
and we have compared our results with existing ones in the literature. For small values of v, 




V 



FIG. 6: Plot of — A c /A as a function of v, for small values of the latter. We have superimposed to 
the data the best fit we could find, given by eq. ((2Tj) with a = 6.3, b = 9.0. 



see figure 6, we find a very good linear fit to the quantity vlog(— A c /A), which means that, 
in that region 

We have compared these results to the estimate given in Ref. Q], i.e. A ~ ee -1 / £ , where 
e = v 2 . This estimate was made for very small values of v and assuming the constraint 
Xv 2 = 1. Taking these facts into account, we find remarkable similarity between that result 
and the numbers we obtain as, in our language, that estimate reads A c /A ~ e^ 1 ^ 2 , to be 
compared to eq. (^1]) . However, as mentioned above, numerical difficulties prevent us from 
exploring the region of very small v and, therefore, we cannot check the validity of this 
equation in that region. 

Having already mentioned the numerical difficulties encountered when exploring the small 
v region, which corresponds to the weak gravity limit, we turned to explore large values for 
v. As it happens, there is an upper bound on the value of A c if we want to obtain regular 
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solutions that trap gravity. This is given by the quantity m\jv 2 vanishing, or, in other words, 
by 

1 + ^-0. (22) 

We show this graphically in figure 7, where we plot — A c / (Av 4 ) as a function of v, indicating 
the points that were calculated numerically. We can clearly see that we are approaching 
the super heavy limit, given by eq. (}2"2"|) . and we can calculate how the different relevant 
quantities in our problem approach it. For example, the asymptotic value of the scalar field, 
/oo, using eq. (fT5j) . will be given by 

f£ = yf v , (23) 

which means that, in the limit where gravity is strongest, the field will be maximally dis- 
placed from the minimum of its potential, at a value given by 



V(f£) = X 1 - ■ (24) 



i 



The other two asymptotic quantities will be given by 

i— v 2 



7Mf 



(25) 



L s h _ V7 1 



oo 



2 V\v 

Note that, in terms of adimensional quantities, = y/Ev/(7Mg), and L s ^ = \fl /2. 

Finally, we should comment on the phenomenological applications of our results. One 
of the motivations to study these kind of set-ups is to try to generate a large hierarchy 
between the higher dimensional and the D = 4 Planck masses in order to explain the so- 
called hierarchy problem of traditional grand unified theories. Here these two scales are 
related by the following equation 

Planck = 5^ / dxM 2 (x)L(x) , (26) 

where, as explained after eq. (11), m(x) = M'(x)/M(x). We have checked that, in the 
solutions we have found, the hierarchy is never larger than a few orders of magnitude (order 
100 for v — 1, increasing to order 1000 for v = 0.7), and obtaining a D = 6 scale of the 
order of the TeV proves numerically unachievable. It would correspond to a region with very 



small v, as already pointed out in Ref 



which we cannot access given that, as mentioned 



before, our problem turns into a two-scale one, and our numerical tools become insufficient. 
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FIG. 7: Plot of -A c /(Av 4 ) as a function of v. 



IV. DISCUSSION AND CONCLUSIONS 



There are in 
six dimensions 



the 



iterature several confining gravity solutions that involve scalar vortices in 
3,0]- In particular, complete numerical solutions have been presented 



in two cases |l9t |20| . It is interesting to compare our work with these solutions. 

In the first one, Ref. Q], Giovannini et al. studied a general gravity trapping abelian 
vortex. The presence of the gauge field is crucial and determines the properties of the 
solution. It is described by three parameters, A, v, and a new one, a, that fixes the gauge 
coupling. The authors find that the condition of gravity trapping translates into a fine tuning 
among these three parameters, fixing a critical surface in this space. The corresponding 
metric, that is a solution to the field equations, exhibits both Minkowski and angular warp 
factors that decrease when r increases. This has to be compared to our case, where only the 
Minkowski factor decreases since the solution is cigar-like. This fine tuning is a generic fact 
in this kind of models, and appears when we try to connect the required behaviours for the 
metric: a regular one at the origin and a confining one for large values of r. The presence 
of the gauge field does not change this fact. 

In the second one, Ref. j^J, a cigar-like solution based on a BPS scalar vortex was studied. 
The potential was derived from a Supergravity (SUGRA) like theory and its structure was 
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almost fixed. In particular, there was no room for a cosmological constant. Nevertheless, 
this SUGRA inspired potential is negative near the local minimum, acting like a negative 
cosmological constant for values of the scalar field close to it. The choice of the potential was 
dictated by simplicity, since for that specific form the field equations are first order (for more 



details see Ref. 



2l|). Therefore numerical solutions are easily obtained just by integrating 



these first order equations together with regular boundary conditions at the origin, which 
happen to be gravity trapping. The SUGRA like structure induces the connection of the 
two (origin and large r) required metric behaviours in the solution. Notice that there is an 
implicit fine tuning in the model in the choice of the potential. If we slightly change the 
coefficients of the terms appearing in the potential, the new model will almost certainly not 
admit a confining solution. This is related to the stability of the solution against radiative 
corrections, which is an issue is beyond the scope of our work. 

To conclude, in this paper we have analyzed the global string in a six- dimensional space 
time with a negative bulk cosmological constant, A. We have presented numerical evidence 
of the existence of solutions that confine gravity. For every value of v, the Higgs vev, there 
is unique value of A that provides a regular solution. This critical cosmological constant is 
bounded by — V(0) < A c < and approaches its lowest value in the strong gravity limit. 
On the other hand, it is difficult to get a hierarchy between Mq and Mpi anc k, at least in the 
range we have been able to explore numerically. 
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